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ABSTRACT 

We report the Suzaku detection of the earliest X-ray eclipse seen in the recurrent nova U Scorpii 
2010. A target-of-opportunity observation 15 days after the outburst found a 27±5% dimming in 
the 0.2-1.0 keV energy band at the predicted center of an eclipse. In comparison with the X-ray 
eclipse depths seen at two later epochs by X MM- Newton, the source region shrank by about 10- 
20% between days 15 and 35 after the outburst. The X-ray eclipses appear to be deeper than or 
similar to contemporaneous optical eclipses, suggesting the X-ray and optical source region extents 
are comparable on day 15. We raise the possibility of the energy dependency in the photon escape 
regions, and that this would be a result of the supersoft X-ray opacity being higher than the Thomson 
scattering optical opacity at the photosphere due to bound-free transitions in abundant metals that 
are not fully ionized. Assuming a spherically symmetric explosion model, we constrain the mass-loss 
rate as a function of time. For a ratio of actual to Thomson opacity of 10-100 in supersoft X-rays, we 
find a total ejecta mass of about 10 -7 -10 -6 M . 

Subject headings: stars: novae, cataclysmic variables — stars: individual (U Scorpii) — X-rays: stars 



1. INTRODUCTION 

A nova explosion occurs in an accreting binary system 
comprising a white dwarf and a red dwarf companion. 
When the amount of accreted material reaches a criti- 
cal mass, hydrogen fusion is trig gered by a thermonu - 
clear runaway on the white dwarf (jStarrfield et al.l feOOS). 
An event is mainly characterized by the development of 
photospheric emission powered by nuclear burning af- 
ter the explosion. This is first in the optical, where a 
sudden increase in brightness is the result of radiative 
transfer through the optically-thick ejecta. As the ejecta 
expand and become less opaque, the dominant emission 
shifts toward higher energies, eventually becoming super- 
soft X-rays. The spectral hardening occurs because the 
mass outflow diminishes with time and consequently the 
pseudo phot osphere is formed deeper in the expanding 
ejecta ( e .g.. iBathl 11978ft. For detaile d reviews see e.g., 
IWarnerl 1)2003 ) : IBode fc EvansI (poM ). 

Our understanding of nova evolution to date has al- 
most entirely resulted from photometric and spectro- 
scopic studies. While all novae are thought to be binaries, 
and some fraction will inevitably be eclipsing, the uti- 
lization of eclipses during an explosion as a powerful and 
direct probe of the emitting geometry has not yet been 
fully realized. Observing such a phenomenon is difficult 
for most novae because they spend many thousands of 
years between outbursts and are therefore observed only 



dtakei@head.cfa.harvard.edu 

1 Smithsonian Astrophysical Observatory, 60 Garden Street, 
Cambridge, MA 02138, USA 

2 Japan Aerospace Exploration Agency, Institute of Space and 
Astronautical Science, 3-1-1 Yoshino-dai, Chuo-ku, Sagamihara, 
Kanagawa 252-5210, Japan 

3 European Space Agency, XMM-Newton Observatory SOC, 
SRE-OAX, Apartado 78, 28691 Villanueva de la Canada, 
Madrid, Spain 

4 Department of Physics and Astronomy, University of Leices- 
ter Leicester, LEI 7RH, UK 

° School of Earth and Space Exploration, Arizona State Uni- 
versity, Tempe, AZ 85287-1404, USA 

6 Department of Physics, Rikkyo University, 3-34-1 Nishi- 
Ikebukuro, Toshima, Tokyo 171-8501, Japan 



once, with no prior information on the nature of the pro- 
genitor. However several novae have been observed to ex- 
plode more than once in a human lifetime and are called 
recurrent novae. Since they can be studied in detail prior 
to an eruption, they are invaluable for understanding the 
nature of novae and cataclysmic variables. 

In this letter, we investigate the development of the 
supersoft source (SSS) region for the eclipsing recurrent 
nova U Scorpii 2010. Using the Suzaku satellite we have 
detected a dimming 15 days after the outburst, corre- 
sponding to the earliest post-outburst X-ray eclipse ever 
reported for this object. Through comparison with a 
30-50% dimming in X-ray eclipse s at two subseq uent 
epochs observed by XMM-Newton ()Ness et alJl2012[ ). we 
find that the X-ray source shrank with time. For the 
first time we are able to use direct geometric information 
provided by supersoft X-ray eclipses to investigate the 
evolution of the radiatively-driven flow in the early stage 
of an explosion. 

2. TARGET (U SCORPII) 

U Scorpii is an eclipsing binary undergoing nova explo- 
sions recurrently about every 10 years. Novae have been 
observed ten ti mes, the most re cent of which were in 
1999 and 2010 ()Schaeferll2010al) . IStarrfield et "all JTH1) 
noted that the short time between outbursts points to 
a mass-gaining white dwarf close to the Chandrasekhar 
limit. The system has an inclination —80 deg, and the 
radii of the companion star and the bi nary orbit are 
2.7 and 6.9 i? Q , respectively, accord ing to lHachisu et al ' 
(2000a), and are 2.1 and 6.5 Rq bv lThoroughgood et al 
(|2001l) . The orbital p eriod before the 1999 out - 
burst was 1.2305521 d (jSchaefer fc Ringwaldl 119951) . 
while the period from 2001 to 2009 was 1.23054695 d 
(|Schaefer et al.ll2010cf). indicating that it ch anged due to 
the 1999 eruption (iMatsumoto et al.1 12003 . We adopt 
the iSchaefer et al.l (|2010cD ephemeris with an origin of 
HJD 2451234.539. An optical eclipse covers -18% of an 
orbital ph ase in the quiescent state, which corresponds 
to 0.22 d (jSchaefer et al.ll2010cD . 
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Fig. 1. — Development of (a) optical and (b) X-ray brightness 
in the nova outburst of U Scorpii 2010. The origin of the abscissa 
is MJD 55224.4385 when the nova was discovered dSchaefer et al.l 
I2010al ). The times of the Suzaku and XMM-Newton observa- 
tions are shown by vertical dashed lines in green and blue, respec- 
tively, (a) Optical magnitudes are from the American Associa- 
tion of Variable Star Observers (AAVSO), Variable Star Observer s 
League in Japan (VSOLJ), and table 2 in Schacfcr et al. (2011). 
(b) Ba ckground-subtracted count rates of the Swift (Schacfcr et al. 
2010b]) in the 0.3-1.0 keV energy band. 

The tenth observed eruption was discovered on 2010 
January 28.4385 UT (MJD 55224.4385) by B. G. Harris 
(jSchaefer et al.ll2010al ). We define the epoch of the dis- 
covery as the origin of time. A worldwide collaboration 
was organized, and subsequent studies wer e conducted by 
ground- and space-ba sed telescopes (e.g.. IQsborne et al.l 
l2010t fSchaefer 2010jj). The development of optical and 
X-ray brightness is shown in Figure [TJ Optical eclipses 
imply the source region was 4.1, 3.4, and 2.2 i?^ on days 
15-26 , 26-41, and 41-67, respectively (Schacfcr et al. 
120111 ). The eclipse shapes were consistent with a spher- 
ical source until day 26, with apparently more disk-like 
morphology at later times. A 30-50% dimming in X-rays 
during expected times of e clipse were foun d by XMM- 
Newton on days 23 and 35 (jNess et al.ll2012j) . suggesting 
the size of the X-ray source was comparable to the optical 
and orbital sizes. The X- ray light curve on day 23 ex- 
hibited oscillations, which iNess et al.l (|2012l ) interpreted 
in terms of a reforming accretion disk. 

3. OBSERVATIONS AND REDUCTION 

We performed target-of-opportunity observations of 
U Scorpii on 2010 February 6, 9, and 12 (9, 12, and 15 
days after the outburst, r espectively) with th e Suzaku 
X-ray satellite. Suzaku ( Mits uda et al.l l2007f) has an 
X-ray Imaging Spectrometer (XI S: iKovama et al.l|2007f) 
and a Hard X-ray De tector (HXD: lTakahashi et alj|2007t 
iKokubun et al.ll2007f ). We concentrate on the XIS data 
on day 15 that show a clear eclipse. The day 9 data 
are poor in statistical quality, and the light curve on 
day 12 exhibited substantial stochastic variability that 
prevented unambiguous measurement of the eclipse sig- 
nature. 

The XIS is equipped with four X-ray CCDs at 
the foci of four co-al igned X-ray telescope modules 
(jSerlemitsos et al.ll2007f ). Three of them (XIS0, 2, and 3) 
are front-illuminated CCDs sensitive in the 0.4-12 keV 
energy band, while the remaining one (XIS1) is a back- 
illuminated CCD sensitive in the range 0.2-12 keV. XIS2 
and a part of XIS0 are not functional and their data were 
excluded. The XIS was operated in the normal clocking 
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Fig. 2. — (a) XIS 0.2-1.0 keV count rates in eclipse (blue) and 
at other times (black), together with the background (red). The 
time of the optical eclipse is illustrated by the gray shaded region 
and the black dashed-and-dotted line. The best-fit quadratic model 
applied to the data outside eclipse is represented by a black dashed 
line. The yellow region shows the ±90% range of the black data 
points. The best-fit eclipse models are shown color-coded by solid 
lines; the quadratic component is illustrated by the dashed lines. 
The vertical dashed-and-dotted line in magenta and cyan indicate 
the Gaussian and transit center, respectively, (b) The hardness 
ratios (HRs) defined by (H— S)/(H+S), where H and S are rates 
in the 0. 5-1.0 keV and 0.2-0.5 keV energy bands, re spectively, (c) 
The V-band magnitudes from Schacfcr et al. ( 2 0111 '). 

mode with an 8 s frame time. 

Data were processed with pipeline version 2.4.12.27. 
Events were removed during South Atlantic anomaly pas- 
sages, when night-earth elevation angles were below 5°, 
and day-earth elevation angles were below 20°. The net 
exposure time is 27 ks. For data reduction, we used 
the HEASoft package version 6.10 and the calibration 
database version xis20090925/xrt20080709. 

4. ANALYSIS 

X-ray light curves and spectra were constructed by tak- 
ing source events accumulated from a circular region of 
130 pixels radius (2 f 3) adaptively-chosen from XIS im- 
ages to maximize the signal-to-noise ratio. Barycentric 
correction of photon arrival times was applied. Events 
taken with the three XIS were merged, and background 
was estimated from an annular region with inner and 
outer radii of 180 and 250 pixels (3fl-4f2), respectively. 
The photon pile-up fraction was <0.1%. Simple spec- 
tral fitting of the extracted photon events on day 15 
yielded similar model parameters (blackbody tempera- 
ture ~ 24 eV; Brems strahlung tempera ture ~ 0.7 keV) 
to those obtained by INess et ail (|2012D for the day 23 
XMM-Newton spectrum, indicating only slow spectral 
evolution through this period. The SSS component dom- 
inates the X-ray spectrum, in which the Bremsstrahlung 
component contributes at most only ~1% of the total 
flux and/or count rates in the 0.2-1.0 keV energy band. 

Effective areas were computed by incorporating the 
history of the satellite aspect solution using a M onte 
Carlo simulation tool fxissim; llshisaki et al.1120071 ) that 
also accounts for an aspect-related vignetting of count 
rates by 10-20% due to the loss of a gyro system on 
Suzaku. Finally, all count rates were normalized to values 
corresponding to extraction from a circular aperture of 3' 
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radius. The resulting background-subtracted light curve 
and hardness ratios in the 0.2-1.0 keV e nergy band, to- 
gether with the F-band magnitudes from ISchaefer et al.1 
poll , are shown in Figure [2] 

A shallow dip in the X-ray light curve is coincident with 
the predicted optical eclipse on day 15 (Figured). Com- 
parison of the distribution of counts during the optical 
eclipse with those outside of eclipse yields a null hypoth- 
esis probability that the X-ray dimming is explained by 
a random source fluctuation on top of the global trend 
of <0.01. No significant correlation was found between 
the flux variation and the spectral hardness (Figure Gp), 
suggesting that the variation of intrinsic colors does not 
affect the dimming. We thus conclude that the X-ray 
dimming is due to the binary eclipse. 

The X-ray eclipse was measured by fitting the entire 
light curve using a quadratic function plus a Gaussian 
component to represent the source emission and approx- 
imate the eclipse dimming, respectively (Figure^). The 
model yielded an eclipse depth of 27±5% at a center of 
15.006±0.012 days after the outburst with a duration of 
0.17±0.04 d (0.13±0.03 orbital phase) in the full width at 
half maximum (la errors). The eclipse center was shifted 
by ^0.03 d in compar ison with the optical ephemeris of 
ISchaefer et al.1 (|2010cD . At face value this is a significant 
offset. We n ote, however, that the optical eclipse minima 
observed bv lSchaefer et al.1 ([201 ID show jitter of a similar 
magnitude. This jitter is likely due to small fluctuating 
brightness inhomogeneities and departures from perfect 
sphericity in the emitting region. 

5. DISCUSSION 
5.1. Source Emitting Region 

The X-ray eclipse provides a means for constraining the 
geometry of the emitting region. Assuming a spherically- 
symmetric source with no limb- darkening and the sys- 
tem parameters of lHachisu et al.l (|2000a|) . the 27±5% 
dimming on day 15 implies that the X-ray source had 
roughly four times the area and twice the radius of 
the companion star, corresponding to a source radius of 
5.1±0.6 Rq. This is consistent with a fit to the light 
curve using a quadratic fun ction and analytic transit 
model (IMandel fc Agolll2002T) assuming the same geom- 
etry, for which we obtained a radius 4.5±0.2 Rq and 
center epoch 14.988±0.007 days post-outburst; see Fig- 
ure [2^,. By ad opting the slightly m o re com pact system 
parameters of iThoroughgood et al.l (|2001| ). radius esti- 
mates are smaller by a factor of ~0.8. This emitting 
region within the residual outflow is large in comparison 
to the unde rlying white dwar f and was referred to as a 
"corona" bv lNess et all (|2012| ). 

The development of the source region can be studied 
by comparing the above radii with the XMM-Newton re - 
sults taken at two different epochs (jNess et al.l I2012T ) . 
The X-ray light curve on day 23 exhibited a dimming 
of up to ^50% with ^20% oscillations during an eclipse. 
The latter coul d be caused by ab sorption in a reforming 
accretion disk (N ess et al.l I2012T ) . or by the geometrical 
eclipse of a photometrically varying photon escape region 
caused by inhomogeneities and/o r density instabilities in 
the radiatively-driven flow (e.g. , iShaviv fc Dotanl 120101 : 
lShavivll2005l: lOwocki et al.lll988D . In contrast, the X-ray 
light curve on day 35 showed a clear eclipse with a ~40% 



dimming, corresponding to a source radius of ~4.2 
indicating that the X-ray source shrunk by about 10-20% 
between 15 and 35 days after the outburst. This is the 
first quantitative X-ray measurement of the shrinking of 
the source radius of an emerging SSS in a nova outburst. 

The standard nova paradigm posits that the SSS emis- 
sion emerges as ejecta expansion and decline in the 
radiatively-driven flow allows photon escape at succes- 
sively smaller radial distance from a white dwarf surface 
until supersoft X-ray temperatures are reached. Once 
this begins and the outer ejecta are thin, the source ra- 
dius is of the order of a few solar radii or smaller and 
the atmosphere is po ssibly puffed-up and porous (e.g., 
IShaviv fc DotanlfeOlOft . This allows it to be much larger 
than the Eddington luminosi ty and X-ray temperature 
would otherwise suggest ('e.g.. lShavivll2005l) . A few thou- 
sand km/s outflow traverses the source in only ^100 s; 
at this point, any further change in radius probes the in- 
stantaneous mass-loss rate in the outflow. We investigate 
this below, first assuming the simplest ejecta model, and 
then examining the implications of an energy- dependent 
opacity by comparing source radii in optical and X-rays 
on day 15. 

5.2. Simple Ejecta Model 

Assuming a spherically-symmetric flow from the white 
dwarf surface, the X-ray photon escape radius R x is re- 
lated to the optical depth by 



Ron 



npdr, 



(1) 



where R ou t is the outermost radius of the ejecta, k is the 
ejecta opacity, and p is the density as a function of radial 
distance r. We define the X-ray photon escape layer as a 
point at which r is approximately unity. Assuming that 
the ejected gases are accelerated within a small radial dis- 
tance and that the terminal velocity Woo has not changed 
appreciably since the eruption (i.e., the time scale of the 
terminal velocity decline is much larger than the time for 
the ejecta to reach R ou t), Rout ~ Vcot. We adopt Voo = 
3000 km s _1 from the estimate of lYamanaka et al~l (|2010T ) 
based on Ha line profiles, and k = 0.4 cm 2 g _1 as the 
Thomson scattering opacity in a fully-ionized hydrogen- 
dominated gas. The radial density profile is given by 



1 M 



(2) 



where M is the mass-loss rate. We next assume a time- 
depe ndent mass-loss rate from the white dwarf given by 
(e.g.. IBode fc Evans|[200l 



M = M Q (t /t) p , 



(3) 



where Mq is an initial value normalized at time t = to = 
1 s and p is a measure of the speed of its decline. Using 
these relations, we fitted the observed radii of the X-ray 
source regions at three epochs (days 15 [Suzaku], 23 and 
35 [XMM-Newton]) with two free parame ters Mp and p. 
We use d the Gaussian estimate with the lHachisu et al.1 
(2000aj) parameters on day 15 in order to compare di- 
rectly the eclipse depths with the other epochs. We 
further adopted 10% uncertainties on the source radii 
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Fig. 3. — Top: The best-fit source radius model vs. time (solid) 
and la confidence (dashed) for radii inferred from Suzaku and 
XMM-Newton eclipses. Models for day 23 eclipse depths of 30% 
(red, case 1) and 50% (green, case 2) are shown separately. The or- 
ange region and the dashed line indicate the sizes of the companion 
star and the orbital separation, respectively. Bottom: Total ejected 
mass vs. time corresponding to the upper panel models. The ver- 
tical dashed-and-dotted line shows the end of the SSS phase. 

derived by XMM -Newton based on the light curves in 
INess et all (|20H . The source radius on day 23 depends 
on the interpretation of the light curve oscillations, thus 
both cases of 30 and 50% eclipse dimming (cases 1 and 
2, respectively) were treated separately. The best-fit re- 
sults and la confidence regions are shown in Figure [31 
where M = lxl0~ 3 (6x 10~ 5 -3x 10~ 2 ) M yr" 1 and p 
= 0.23±0.22 in case 1; M = 5xl0~ 3 (3x 10~ 4 -2x lO" 1 ) 
M Q yr" 1 and p = 0.33±0.25 in case 2. 

Integration over the time-dependent mass-loss rate 
yields a rough estimate of the total ejecta mass M ej - as 
a function of elapsed time and is shown in Figure [3] As- 
suming that the SSS faded and the mass outflow was 
terminated on day 40 (Figure Qi>) , the total ejecta mass 
is about 5-8 xlO -6 Mq for the simple model. This is con- 
sistent with the changing orbital period from long-term 
optical monitoring (4.3±6.7x 10~ 6 M Q : lSchaeferil2lnTl ). 
but slightly higher th an theoretical estimat es for U Scor- 
pii (e.g., 4xlQ- 7 : IStarrfield et al.1119881 : 2xlO- 6 M : 
IHachisu et al.ll2000bD . We here note that mass-related 
values al so become smaller by the factor of ^0.8 by 
adopting iThoroughgood et al.l ([2001) parameters. 

5.3. Optical vs. X-ray Eclipses 

ISchaefer et al.l (|2011D found that the eclipse shapes in 
optical were consistent with a spherical source with a ra- 
dius of 4.1 Rq on days 15-26. Since we have assumed 
spherically-s ymmetric source reg ions whereas the opti- 
cal eclipses (fechaefer et al.ll201lD and the day 23 XMM- 
Newton eclipses ()Ness et al.H2012f) were accompanied by 
more complicated structure, we cannot conclude with 
certainty that there is a significant difference between 
the sizes of the optical and X-ray photon escape re- 
gions, though at face value the different inferred radii 
suggest this on day 15 with the system parameters of 
IHachisu et all (|2000aD . 

A remaining puzzle is the origin of the optical plateau 
started from ~10-15 days post-outburst. This opti- 
cal flux clearly exceeds the extrapolated X-ray black- 
body with the Eddington luminosity of a Chandrasekhar 
white dwarf, though a blackbody is grossly inadequate 
to describe the true spectrum and the magnitude of 



any op tical excess remains uncertain. IHachisu et al.l 
(2000a) argued the optical flux included reprocessing 
photospheric emission on the surface of the reappear- 
ing accretion disk. Currently, no clear evidence is found 
for this because it is unknown how early the accretion 
process resumes after the explosion, though recent stud- 
ies imply that the non-accreti on case is more likely, at 
least on day 15. Simulations of iDrake fc Orlandol (|2010f ) 
found that the accretion disk was c omple tely destroyed 
by th e blast. IMason et all (|2012D and IWorters et all 
(2010) suggest resumption of accretion based on spec- 
troscopy and optical flickering, respectively, around day 
8, t hough the short term variability disappea red on day 
15 (jMunari et al.|[20Tol) . ISchaefer et all (poll found the 
optical so urce spherical un til day 26 based on eclipse 
mapping. Nes s et al.l (|2012|) further discussed reforming 
accretion still on day 23. Alternatively, the reprocessing 
might occur in surrounding ejecta with a substantial op- 
tical depth in supersoft X-rays, though it has to be noted 
that any scenario at this time remains largely specula- 
tive. 

If we assume that the optical flux on day 15 origi- 
nated mostly from the surrounding ejecta, the possible 
difference between the radii of optical (R ~ 4 Rq) and 
X-ray (R x ~ 5 Rq) sources can arise from the energy 
dependence of the photon escape radii that would result 
from different optical and X-ray opacities (n and n x , re- 
spectively). For a plasma at temperatures of 10-100 eV, 
the opacity in the supersoft X-ray range can be signif- 
icantly larger than the Thomson cross-section owing to 
abundant elements (primarily C, N, O, and Ne) not being 
fully ionized. For a low-density plasma, the optical opac- 
ity will instead be close to the Thomson value. Using th e 
PINTofALEQ routine IONABS (jKashvap fc Drakd2000l) . 
we computed the cross-section assuming collisional equi- 
librium at the temperature of 50 eV and found k x ~ 
50 cm 2 g~ x in the 0.2-1.0 keV range, or ^100 times the 
Thomson value. The plasma in the ejecta is in fact likely 
to be photoionized rather than in collisional equilibrium. 
For a constant wind velocity Eqns. Q] and [2] imply for op- 
tical depth unity npr ~ 1, or p ~ 5xl0 12 H atoms cm' 3 
at 4 Rq. T he ionization parameter, £ = L/n e r 2 (e.g., 
iTarter et al1l!969f ). is then of the order of 10 2 for the Ed- 
dington luminosity and indicates the outflow is photon- 
dominated. Nevertheless, our cross-section estimate sug- 
gets that the true supersoft X-ray opacity could be an 
order of magnitude or more larger than that in the opti- 
cal. 

The density profile, and consequently the mass-loss 
rate and total mass loss, derived from the X-ray eclipses 
scale inversely with the gas opacity, so that the total 
mass-loss rate would be 5-8xl0~ 6 {k /k x ) Mq. We 
have argued that the X-ray opacity at the photosphere is 
higher than the Thomson value that characterizes the op- 
tical range, and that this would give rise to the observed 
difference in optical and X-ray source radii. Theoretical 
mass-loss estimates are in the range 4x 10~ 7 - 2 x 10 -6 Mq 
(|Starrfield et al.lll988t IHachisu et aLll2000bD . indicating 
the true X-ray opacity might b e higher than the optica l 
one by a factor of about 2-20. IDrake fc Orlandol 1)20101 ) 
found the initial explosion threw off probably no more 
than 10" 7 Mq using non-spherical hydrodynamic mod- 

7 http: / /hea- www. harvard.edu/PINTofALE/ 
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els and early X-ray luminosity constraints. Estimates of 
the subsequent mass loss then imply that most of the 
mass was lost during la ter evolution, confirming: the oret- 
ical expectations (e.g.. iGallagher fc Starrfieldlll978t ). 

We thank Suzaku for undertaking our target-of- 
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